Self-propulsion of V-shape micro-robot 



CO 



By V. A. Vladimirov 

Dept of Mathematics, University of York, Heslington, York, YO10 5DD, UK 

(Received Sept 12th 2012) 

In this paper we study the self-propulsion of a symmetric l/-shape micro-robot (or 
T^-robot) which consists of three spheres connected by two arms with an angle between 
them; the arms' lengths and the angle are changing periodically. Using an asymptotic 
procedure containing two-timing method and a distinguished limit, we obtain analytic 
expressions for the self-propulsion velocity and Lighthill's efficiency. The calculations 
show that a version of l^-robot, aligned perpendicularly to the direction of self-swimming, 
is both the fastest one and the most efficient one. We have also shown that such F-robot 
is faster and more efficient than a linear three-sphere micro-robot. At the same time 
the maximal self-propulsion velocity of V-robots is significantly smaller than that of 
comparable microorganisms. 
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1. Introduction and formulation of problem 

The studies of micro-robots represent a flourishing modern research topic which strives 
to create a fundamental base for modern applications in medicine and technology, see e.g. 



Purcell (1977) Becker et.al. (2003) Najafi fc Golestanian (2004) Dreyfus et.al. (2005) 



Earl et.al. (2007)[|Chang et.al (2007)[|Alouges et.al. (2008)[|Gilbert at.al (2010)[|Golestanian fe Ajdari (2008) 



Golestanian fc Ajdari (2009)||Alexander et.al. (2009)[|Leoni et.al. (2009)[|Lauga (2011) 
The simplicity of both time-dependence and geometry of micro-robots represents their 
major advantage in contrast with extreme complexity of self-swimming microorganisms 



e.g. |Pedley fc Kessler (1987)||Vladimirov et.al. (2004)||Pedley (2009)[|PoIin et.al. (2009) 
This advantage allows to describe the motion of micro-robots in a greater depth. In this 
paper, we study the self-propulsion of a symmetric y-shape three-sphere micro-robot 
(for brevity we call it V- robot), see the figure. We employ two-timing method and distin- 
guished limit arguments, which lead to a simple and rigorous analytical procedure. The 
self-propulsion velocity and Lighthill's swimming efficiency are calculated analytically. It 
appears that V-robot can swim in both directions of x-axis, and V-robot, aligned per- 
pendicularly to the direction of self-swimming, is the most efficient one. In addition, such 
y-robot is faster and more efficient than a linear three-sphere micro-robot, the most stud- 
ied one. y-robot has already been studied numerically by Earl et.al. (2007)| but never 
analytically. Our approach is technically different from all previous methods employed 



in the studies of micro-robots (except Vladimirov (2012b)). The possibility to derive an 
explicit formula (|2.10[) for a y-shape three-sphere micro-robot shows its strength and 
analytical simplicity. The used version of the two-timing method has been developed in 



Vladimirov (2005) Vladimirov (2008) Vladimirov (2012a) 



We consider a symmetric ^-shape micro-robot (V-robot) consisting of three rigid 
spheres of radii R v , v = 0,1,2 (i£j = R 2 ) connected by two arms of equal length I. 
The angle between the arms is 2tp. In plane cartesian coordinates the centers of the 
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spheres and the distances between them are 

We accept the notations (see the figure): 

x (o) = ( a; (o) )2/ (o) ) = (X) o) ) *W = (s^yC 1 )) = (a, y) = (* + £,!,), (1.1) 

r( 01 ^(£,y), r<° 2 > = r< 21 ) = (2 V , 0); 

| p (oi)| = | p (02)| =i] | r (2i) | = 2 y = 2Zsin ¥ >, (£,y) = I (cos ip, simp); 

n (oi) _ (cosy;, sin yj), n*- 02 - 1 = (cos — s'mip), nP 1 " 1 = (0,1) 

where are unit vectors along r^ v \ We use latin subscripts (i, k = 1, 2) for cartesian 

components of vectors and tensors, subscript a = 1,2,3 - for generalised coordinates, 
and subscripts (or superscripts) n, k< = 0, 1, 2 to identify the spheres. F-robot represents 
a mechanical system described by three scalar parameters (generalised coordinates) 

q=(qi,q2,q 3 ) = (X,l,cp) (1.2) 

y-robot moves due to the prescribed motion of the arms 

1 = L + eI(t), Lp — <f> + £<J5(t); t = cut, uj = const, e = const (1-3) 

where the functions Z(r) and £>(t) are 27r-periodic with zero average values, while L and 
$ are constants. The spheres experience external friction forces = (F[ v \ F^); how- 
ever, the arms are chosen to be so thin in comparison with any i?„, that their interaction 
with a fluid is negligible. 

The considered problem contains three characteristic lengths: the length of arms £ c har, 
the radius of spheres i? c harj and amplitude of the arm's oscillations a c h a r- The character- 
istic time-scale is T char = 1/cj and the characteristic force is F char . We have chosen 

ichar = L, i? c har = (Ro + 2i?l)/3, &char — 

sL, 

^char 

= 67T?7i? c har£char/T , c har(l-4) 

where rj is viscosity of a fluid. Two small parameters are 

£ = flchar/L < 1, 6= 3i? ch ar/(2L) < 1 (1.5) 
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The dimensionless (asteriated) variables are chosen as 

X = Lx* , Ri = R c h ar R*, t = T c h ar t*, fi = -Fckar/*! (1.6) 

Below we use only dimensionless variables, however all asterisks are omitted. 

Generalized coordinates q = q(t) (|1.2[) determine the motion of IZ-robot. It can be 
described by the Lagrangian function £ = C(q, q t ), which includes the constraints (|1.3p 

C(q, q t )=tC + f(l-l- el) + g(cp - $ - e<p) (1.7) 

where subscript t stands for d/dt, f and g are Lagrangian multipliers, and K. is kinetic 
energy of a robot; notice that in dimensionless variables L = 1. The Lagrange equations 
are 

dtdq a t oq a dq ° 

where Q = (Qi, Q2, Q3) is the generalized external force, exerted by a fluid on T^-robot. 

The fluid flow past F-robot is described by the Stokes equations. In consistent approx- 
imation masses of spheres and arms are negligible, then /C = 0. Hence (|1.7|) . (|1.8p produce 
the system of equations 

Qi=0, Q 2 + f = 0, Qs+9 = (1.9) 
The calculations of Q a (|1.8p with the use of (jl.l[) yield 

Qi=F 1 (0) +F+, Q 2 = F + cos tp + F~ simp, (1.10) 
Q3 = — IF + sin <p + IF~ cos ip; 

The substitution of (| 1 . 10[) into (|1.9p and simple transformations yield 

Fl 0) + F+ = 0, F+ = A, F ~ = a; (1.11) 
A = —/cos ip + {g/l) simp, a = — {g/l) cost/3 — / simp, 

The explicit expressions for F^ arc 

F M ~ -R v x { t ) + 8 R^RuS^x^ (1.12) 

S ~^fc - 1^)1^ + ™i ™fc J 

Each force i* 1 ^ represents the first approximation for the Stokes friction force exerted on 
a sphere moving in the flow generated by the other two spheres. To construct (|1.12p we 



use a classical explicit formula for a fluid velocity past a moving sphere, see Lamb (1932) 



Landau fc Lifshitz (1959)] |Moffatt (1996)] The substitution of (fI~T2"ll into (fTTTTj) yields: 



Rox^ + Ri(x$ + x$) -SA = (1.13) 
Ri(x$ + x$) - SB = A 
R 1 {x$-x%>)-8C = a 
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where 



A — RqRi 



q(01) 



if. 



°lk x kt °lk 

e/i\iy/i 2 ) 



(01) 
k 



o(02) W (0) 
D lk ) x kt 



! c(21) ^(1) 



q(02) 



y(l + m 



D lk 



with summation convention over fc = 1,2. Functions B and C will not affect the calcula- 
tion of the self-propulsion velocity with precision linear in 5, therefore it is sufficient to 
write B = 0(1) and C = 0(1). The equations (fLl3|) . (fO jt . (where definitions (fTTj) should 
be taken into account) represent a system of five equations for five unknown functions of 
time: X, £, y, A, and a. 



2. Two-timing method and asymptotic procedure 

2.1. Functions and notations 
The following dimensionless notations and definitions are in use: 

(i) s and r denote slow time and fast time; subscripts r and s stand for related partial 
derivatives. 

(ii) A dimensionless function, say G = G(s,r), belongs to class I if G = 0(1) and 
all partial s-, and T-derivatives of G (required for our consideration) are also 0(1). In 
this paper all functions belong to class X, while all small parameters appear as explicit 
multipliers. 

(iii) We consider only periodic in t functions {G £ V : G(s, r) = G(s, t + 27r)}, where 
s-dependence is not specified. Hence all considered below functions belong to Vf]I- 

(iv) For arbitrary G 6 V the averaging operation is 



1 



T0 + 27T 



<G} = — y G( S ,r)dT = G(s), Vr (2.1) 

(v) The tilde-function (or purely oscillating function) represents a special case of V- 
function with zero average (G ) = 0. The bar-function (or mean-function) G = G(s) does 
not depend on r. A unique decomposition G = G + G is valid. 

2.2. Asymptotic procedure and successive approximations 

The introduction of fast time variable r and slow time variable s represents a crucial 
step in our asymptotic procedure. We choose r = t and s = e 2 t. This choice can be 



justified by the same distinguished limit arguments as in Vladimirov (2012a) Here we 
present this choice without proof, however the most important part of this proof (that 
this choice leads to a valid asymptotic procedure) is exposed and exploited below. We 
use the chain rule 

d/dt = d/dr + e 2 d/ds (2.2) 

and then we accept (temporarily) that t and s represent two independent variables. 
Unknown functions X, £, y, A, and a are taken as regular series 

X(t,s) = X (T,s)+eX 1 (T,s)+e 2 X 2 (T,s) + ... (2.3) 

and similar expressions for £, y, A, and a. In (|2.3|) we take Xq = 0, which which express 
the basic property of considered solutions: long distances of self-swimming arc caused by 
small oscillations. After the application of (|2.2|) to (|2.3j) we have 

X t =eX lT +e 2 (X 2T + X 0s ) + O(e 2 ) (2.4) 
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and similar expressions for 

The successive approximations of equations (|1.13[) are: 
Terms of order e° — 1: Ao = and <7o = 0; 
Terms of order e : 

3X lT = -2Ri£ lT , X lT + £ lT = Ai/2, y lT = a±/2 
Its average part gives Ai = and W\ = 0, while its oscillating part leads to: 



2R X 



2R X 



(Zcos$ - (psin$) 



(2.5) 
(2.6) 



2^ 

Xi = (l T cos $ — (p T sin <£) 

3 



(Ti = 2(Z T sin $ + <p T cos $); 
We have used that in dimensionless variables (|1.4p . (|1.6|) 

R + 2Rx = 3, R = 3/(1 + 2p), i?i = i? 2 = 3p/(l + 2p), p = ifc/Jfc; (2.7) 
Terms of order e 2 : Here we consider only the first eqn. ([1.13[) . which can be rewritten as 



(l+e/l 2 )(X t+Xt ) + ^-y t 



RxXt 

y 



(2.8) 



where x = X + £. Since wc consider only linear in 5 precision, then in the second 
approximation we should substitute (into the right hand side of (|2.8jl ) the solutions from 
the first equality in (|2.6[) . Further transformations yield 



^0s — 



SRx 



2MR0 ~ 2Rl) (lG lT ) + 2R a {yxHxr) + ^&K lT ) 



(2.9) 



G=(l + e/n/l, H = ^y/l A , K=l/y 



It is instructive to write that in the original general notations (|1.12[) formulae (|2.9I) . (|2.7[) 
can be presented as 



k=X v=0 [i=£v 

where both tilde- functions in the right hand side are from the first approximation in s. 



2.3. Self-propulsion velocity 

The self-propulsion velocity is defined as Vq = X t — e 2 X 0s . Expressions (|2.9j) . (|2.7[) lead 
to the formula 



V = 2 £ 2 Jt/($,p)(^ T }, 



(2.10) 



3p 



4sin$(l - 6pcos 2 $) 



P 



sin 2 $ 



2(1 + 2p)3 

which represents the main result of the paper. Let us discuss it in detail: 

(i) Vo is proportional to the correlation (l(j) T )] without any restriction of generality we 
consider only (l(f> T ) > 0. 

(ii) Further simplification can be achieved if we accept that the oscillations of arms 
arc harmonic 



Z = cos(t + #i), <p = cos(r + 9 2 ); 2(Z<p T ) = sin( 



(2.11) 
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with constant phases 9i, 62; < 9 < tt. In this case the optimal stroke (providing the 
maximum of Vo) is 9 = tt/2, and max((^ r ) = 1/2, when the self-propulsion velocity is: 

V =e 2 5U(<S>,p) (2.12) 

(iii) Vq (|2.10p , (|2.12[) represents a function of two independent variables $ and p for the 
domain < $ < 7r and < p < 00. Due to symmetry Vo{tt/2 + 4>,p) = Vo{n/2 — </>, p) 
with < </> < 7r/2, the backward-oriented (see the figure) ^-robots with $ = 7r/2 — (f> 
and forward-orientated V- robots with $ = 7r/2 + <p swim with the same velocity. Hence 
it is sufficient to study the domain 

< $ < n/2, < p < 00, 0<6»<7T (2.13) 

where the restriction 9 < it appears due to (l<p T ) > 0. 

(vi) A sharp singularity Vo — > +00 takes place for $ — >■ 0. However this limit does not 
have any physical meaning, since small values of $ correspond to collision and overlapping 
of spheres i?i and i?2- In order to avoid such collision and overlapping, one have to accept 
sin$ > 35/2 > SR\. However, a stronger restriction is required in order to provide the 
validity of approximation for a velocity field in p. 121) . where our basic assumption is: 
the radius of each sphere is much smaller than any distance between them; in particular, 
it means that sin$ ^ SRi. In practice one can take, say, sin$ > 5(5, which for S = 0.1 
gives a 'secure' domain instead of (|2.13j) : 

7r/6<$<7r/2, 0<p<oo, 0<9<it (2.14) 

More generally, the study of motion with small values of $ requires the estimation 
of errors for the used approximation; the best way of doing it is computational, see 



Earl et.al. (2007) We do not consider this complex problem and restrict ourselves to the 
consideration of internal maxima of Vo . 

(v) Vo —5- when Ro — > (or p —> 00), which corresponds to the well-known fact that 
a dumbbell with oscillating arm is not able to swim; also Vo — > when R\ — > (or 
p — > 0), which represents a limiting case of one sphere without any oscillations. 

(vi) A local maximum of Vo for any p = const always takes place at $ = tt/2, when 

U(n/2, p) = 3p(4 + p)/[2(l + 2pf] (2.15) 

Hence, a completely 'open' l^-robot (with 2$ = it) swims in positive direction (see the 
figure) with the maximal speed. Function U(ir/2,p) (|2.19p is increasing linearly for small 
p; U reaches m&xU ~ 0.47 at p ~ 0.27, and then U decreases to zero rapidly and 
monotonically. Hence we can write that 

maxF (<I>, ,0,(9) ~ 0.47e 2 <5 at $ = tt/2, p ~ 0.27, 9 = tt/2 (2.16) 

It is interesting to compare this result with the result for a homogeneous linear three- 



sphere micro-swimmer, where maxl^o — 0.19e 2 5 (see Vladimirov (2012b) I. The compar- 
ison shows that V-robot swims about 2.5 times faster than a linear micro-robot (when 
in both cases the strokes are harmonic and optimal). The fastest swimming of y-robot 
takes place when the spheres R\ and R2 are approximately 4 times smaller than Rq. 

(vii) V- robot can swim in negative direction (see the figure) at $ = tt/4, when the 
function [/(vr/4, p) is non-monotonic and changes its sign: U increases linearly for small 
p, then U reaches m&xU ~ 0.19 at p ~ 0.13, after that decreases such that U = at 
p = 0.43. For p > 0.43 we have U < and min U = f7(1.75) ~ -0.24 which yields 

min y (7r/4,P, tt/2) - -O- 24 e2(S at P ^ 2.10 (2.17) 
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It follows from (|2.10|) that min Vq in all domain (|2.13l) is 

mmV (<S>,p,8) ~ -0.26 e 2 <S at $ = 0.69, p ~ 1.65, = tt/2 (2.18) 

which is close to (|2.17|) . Hence V-robot swims in negative direction if the angle between 
arms is close to 7r/2 (or $ = 7r/4). The maximal speed of this reverse swimming is 
approximately two times slower than the maximal speed in positive direction. 

(viii) It is also of interest that for three equal spheres (p = 1) we have 

maxF o ($,l,0) ~ 0.28e 2 5 at $ = tt/2, 6 = tt/2 
minVo($,l,0) ~ -0.23 e 2 (5 at $ = 0.68, 9 = tt/2 

which shows an essential reduction of the speed in positive direction. 

(ix) The calculations at the boarder $ = 7r/6 of a 'secure' domain (|2.14p show that 
U (7r/6, p) increases linearly for small p, reaches max U ~ 0.13 at p ~ 0.13, then decreases 
such that C/(0.40) ~ and U < for p > 0.40; mint/ = [7(1.65) ~ -0.18. These values 
show that $ = 7r/6 is still well away of a singularity at $ — > 0. 

2.4. Power and efficiency 
The power of V-robot is defined as 

2 

V = f (v) ■ x i U) ( 2 - 19 ) 

where is the force, exerted by the arms on the ^-th sphere. The total force exerted 
on each sphere must be zero, hence + = with the friction force (|1.12p . 
Therefore the main term of (|2.19| can be presented as 

V * £ Ru x<"> ■ * e 2 £ Ru x<& ■ x<& (2.20) 
where we have used (|2~4")l . Then the use of (|2~5|) . (fTTj) . (|2~T |) yields 

15 ~ (lf^) ^ + 2P( ^ + ^ )} (2 - 21) 
and the use of (]2.6[) leads to 

+ iff. + 2p[7 2 sin 2 $ + ^ 2 cos 2 $ + sin(2$)]) (2.22) 



(l + 2p) 2 

For harmonic oscillations (|2.1ip it gives 
— 6e 2 p 



(l + 2p) 2 



W, W = [1 + p(l + sin(2$)cos0)] (2.23) 



t2 



Another expression V s = 3V represents the power, which is required to drag y-robot 
with velocity Vo in the absence of its oscillations when the main approximation for the 
dimensionlcss Stokes friction force is — 3Vo (where the coefficient 3 represents the sum 
of all radii (j2.7jl ). LighthhTs swimming efficiency (see Becker et.al. (2003)) is the ratio 
£ = V s /V. For harmonic oscillations (|2~TT|) expressions ([2TT2]) . (|2T2"3f . (|2~TU|) give 

£ = E^p,e),e^ w f w ^s^e (2.24) 
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The analysis and the computations show that 

maxf ~ 0.89e 2 <5 2 , at $ = tt/2, p ~ 0.15, 9 = tt/2 (2.25) 

One can see that maxf takes place at the same $ = tt/2 and = tt/2, but for ap- 
proximately twice smaller p than for max Vq (|2.16l) : hence the most efficient swimming 
takes place when the spheres R\ and R2 are approximately 7 times smaller than Rq. 
It is also interesting that for a linear three-sphere micro-robot max£ ~ 0.18e 2 i5 2 (see 



Vladimirov (2012b) I, which is five times lower than (|2. 251) , 



3. Discussion 

(i) V-robot has already been studied numerically by Earl et.al. (2007)] but never an- 
alytically. Quantitative comparison between our analytical results and computations of 



Earl et.al. (2007) is difficult, since these authors studied the large amplitudes of arms' 
oscillations and non-harmonic strokes; also the domain of the main parameters, they 
studied, is unclear from the text. At the same time the existence of both forward 
(|2.16p and reverse (|2.18|) swimming provides qualitative agreement of our results with 



Earl et.al. (2007) 



(ii) One can see that the magnitude of velocity in terms of small parameters Vq = 
0(e 2 S) (|2.10|) is the same as the result b y|Golestanian fc Ajdari (2008)||Golestanian fc Ajdari (2009) 



for a linear three-sphere robot and by Vladimirov (2012b) for a linear TV-sphere robot. 
At the same time our choice of slow time s = e*t (|2.2[ ) agrees with classical studies of self- 
propulsion for low Reynolds numbers, see Taylor (1951)| Blake (1971) Childress (1981) 



as well as the geometric studies of Shapere & Wilczek (1989) 



(iii) In our examples, all arms move harmonically (|2.11l) ; it does not provide the maxi- 
mum of Vq ~ ( Hp T ) (j2.10J) . see the relevant discussion in Vladimirov (2012b) The studies 

of non-periodic oscillations represent an interesting additional problem, see Golestanian & Ajdari (2009) 

(iv) In our study we build an asymptotic procedure with two small parameters: e — > 
and 5—^0. Such a setting usually requires the consideration of different asymptotic paths 
on the plane (e, 5) when, say S = 5(e). In our case we can avoid such consideration, since 
small parameters appear (in the main order) as a product e 2 5. 

(v) The mathematical justification of the presented results by the estimation of an error 



in the original equation can be performed similar to Vladimirov (2010) Vladimirov (2011) 



It is also possible to derive the higher approximations of Vq, as it has been done by 
Vladimirov (2010)] Vladimirov (2011) for different cases. The higher approximations can 
be useful for the studies of motion with Vq = (which are always possible for V- robots). 

(vi) In order to compare the velocities of micro-robots and micro-organisms we use 
the dimensional variables, in which maxy o ~ 0A7u!*L*e 2 5; it shows that F-robot can 
move itself with the speed about 10 percent of its own size per second (we have taken 
e = S = 0.2 and uj* = 30s _1 ; the value of w* can be found in Pedley & Kessler (1987) 
Vladimirov et.al. (2004)] |Pedley (2009)] |Polin et.al. (2009)| From these papers we also 



can see that this estimation of V is about 10 times lower than a similar value for 
natural micro-swimmers. 



The author is grateful to Profs. 
Pitchford for useful discussions. 
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